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Abstract 

We show that our procedure of constructing excited-state energy functionals by splitting k-space, 
employed so far to obtain exchange energies of excited-states, is quite general. We do so by applying 
the same method to construct modified Thomas-Fermi kinetic energy functional and its gradient 
expansion up to the second order for the excited-states. We show that the resulting kinetic energy 
functional has the same accuracy for the excited-states as the ground-state functionals do for the 
ground-states. 
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I. INTRODUCTION 



Since the inception of ground-state density functional theory (DFT) 



ASSIES 



5|, efforts 



have been made to extend it to excited-states. Such attempts include the work of Ziegler et 



von Barth [8[, Perdew and Levy [9|, Pathak |10|, T heophilou 



al. Q|, Gunnarsson et al. 

Oliveira, Gross and Kohn Nagy [l^], Sen [15|] and Singh and Deb [16|. However, 



excited-states, akin to its 
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25| only 



a general excited-state density functional theory for individua 
ground-state counterpart, has started taking shape H Q, Q, 
over the past decade or so. 

In density functional theory, energy of a system is expressed as a functional of the density 



of the system. History 26( of writing energy of a system in terms of its density is as 
old as quantum-mechanics itself. In an attempt to simplify the problem of interacting 
electrons, Thomas 27| and Fermi 28[ expressed the kinetic energy of a many electron system 



approximately by employing the expression 



29j for the kinetic energy of the homogeneous 



electron gas (HEG). Similarly, Dirac 30] gave an approximate expression for the exchange 
energy of a many-electron system by employing the corresponding HEG formula. With the 
Hohenberg-Kohn [l] discovery of the one-to-one map between the ground-state density and 
the Hamiltonian of a system, it became clear that the energy of a system can indeed be 
expressed as a functional of its ground-state density; however, the functional is not known 
exactly. In the Kohn-Sham formulation 2| of DFT, the kinetic energy component of the total 
energy is treated highly accurately by writing it in terms of orbitals of an auxiliary system. 
Thus the non-interacting kinetic energy is expressed in terms of the Kohn-Sham orbitals 
\<pi) as (atomic units are used throughout the paper so that we take h — m e — \e\ — 1) 

1 
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(i) 



where /j represent the occupation of i th orbital. For the ground-state, /, = 1 for i th orbital of 
each spin if i < i max where i max is the index of the uppermost occupied orbital and is for all 
the higher orbitals. For an excited-state, the occupation is different from the ground-state; 
for example it could be equal to 1 for i < ii, for i x < i < i 2 and 1 again for i 2 < i < i-3, as 
shown schematically in Fig. 1. On the other hand, the exchange and correlation energies are 
still expressed approximately in terms of the density. Foremost among these approximations 
are the local-density approximation (LDA) and the local spin-density approximation (LSD). 
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In these approximations, the exchange and correlation energies are expressed in terms of 
the density by employing the corresponding expression for the HEG. Thus the LDA for the 
exchange energy is the same as the Dirac expression for it. Over the years, far more accurate 



functionals 3l| for exchange and correlation energies have been constructed by going beyond 
the LDA and including corrections in terms of the gradient of the density. The leading term 
in most of these functionals is the LDA/LSD functional and in the limit of the gradient of 
the density vanishing, the functionals indeed reduce to the latter. 

Given this background, a question that arises naturally in the development of excited- 
state DFT is if it would possible to construct energy functionals for these states with similar 
accuracy as is obtained in the ground-state functionals. In particular it is important to 
develop an LDA functional for the excited-states since that is the foundation on which more 
accurate functionals are built. We have recently constructed an exchange energy functional 
for excited-states within the LDA. This has been done by splitting the k-space in accordance 
to the occupied and unoccupied orbitals of the excited-state, as shown in Fig. 1. In the 
figure, we have some orbitals - the core orbitals - including the lowest energy orbitals that 
are occupied, then some empty orbitals and then some more orbitals - the shell orbitals - 
that are occupied again. The k-space, also shown in the figure, is accordingly split such that 
it is occupied from k = to k = ki, empty from k\ to &2 and then occupied again from &2 
to k 3 . Here k±, k^ and k% are given by the equations 

fc?(r) = 37r 2 p c (r) (2) 

A; 2 3 (r) - fc»(r) = ^ 2 p vac {v) (3) 

kl(r)-k!(r) = 3n 2 p s (r) (4) 

where p c and p s are the electron densities corresponding to the core and the shell orbitals. 
Similarly, p vac is the electron density corresponding to the set of unoccupied orbitals. Thus 

Pc(r) = EjC e (r)| 2 
Pvac(v) = Ej0r° cc W| 2 
P*(r) = £j0f e "(r)r 



3 



Ground state 



-ti- 



an excited state 



(a) 



(b) 



(c) 



FIG. 1: Orbital occupation of electrons (a) and the corresponding /, for each spin drawn contin- 
uously as a function of orbital energy ej (b) for the ground and an excited-state of a finite system. 
The corresponding k-space occupation (c), in the ground and an excited state configuration similar 
to that shown in (a) for a homogeneous electron gas (HEG). 



The total electron density p(r) is given as 

p(r) = p c (r) 



Ps\T) 



(5) 



For detailed derivation of these equations, we refer the reader to the next section. Em- 
ploying the exchange energy functional developed by us, we have been performing accurate 



calculations 
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241 ] of excited-state energies of a variety of systems including the 



band gaps [32| of a wide variety of semiconductors in the recent past. 

Although the results obtained by us with the excited-state exchange energy functional 
are impressive, the question that we have been asking ourselves is if the method employed 
by us - that of splitting the k-space - to construct the functional is general. If the answer 
is in the affirmative, the same method should also lead to reasonably accurate functionals 
for other components, viz. the kinetic and the correlation energies, of the total energy. 
Further, we should be able to build on the LDA to include higher order corrections in terms 
of the gradient of the density. In this paper we address this question in connection with 
the non-interacting kinetic energy of a system of electrons. Our aim in these investigations 
is to explore if a kinetic energy functional for excited states, constructed by splitting the 



k-space, gives similar accuracy for exact kinetic energy of these states as the well known 
Thomas- Fermi or the gradient-expansion approximation (GEA) functionals!], 4] do for the 
ground states. We show in this paper that it does. Thus our present results demonstrate 
the robustness of our procedure of constructing energy functionals for the excited-states of 
a many-electron system. 

We note that besides the__GEA functional, there is another approach to constructing 

for the ground-state, which employs two exact asymptotic 



kinetic-energy functionals 
forms: Thomas-Fermi for the HEG and the von-Weizsacker term 



34j for one-orbital systems. 



If our approach has universality, it should also work for functional such as the Gazquez- 
Robles functional proposed in reference 33(. We show in this paper that it does. 

We start in the next section with a description of the Thomas- Fermi approximation for 
the non-interacting kinetic energy for the ground-state. This approximation is the LDA 
for the kinetic energy. We then discuss the gradient expansion approximation (GEA) for 
the kinetic energy up to the second-order in the density gradient. Results for a few atomic 
systems and the key features of these results are then discussed. This forms the background 
against which the kinetic energy functional for the excited-state is then constructed and 
tested in the section after the next one. We end the paper with some concluding remarks. 

We point out that our aim in this paper is to explore conceptually if our approach yields 
kinetic energy functionals that have accuracy similar to their ground-state counterparts. 
Our work shows that it does. Question arises: Can these functionals be applied to obtain 
average excited-state energies. This possibility is being explored. However, density-based 
functionals cannot be expected to reproduce the exact answer for the kinetic energy, as given 
by equation 1. Therefore the operational utility of the kinetic energy density-functionals for 
excited-states is similar to that of traditional Thomas-Fermi functional or its extensions for 
the ground-state. 



II. LDA AND GEA UP TO THE SECOND ORDER FOR THE NON- 
INTERACTING KINETIC ENERGY OF THE GROUND-STATE 

The basis of the LDA is the homogeneous electron gas for which the kinetic and the 
exchange energies can be expressed in a rather simple form involving the density of the 
system. For the non-interacting kinetic energy we consider a gas of non-interacting electrons 
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that fill the k-space from k = to k = kp because of the Pauli exclusion principle. The 
wavefunction for an electron in a state specified by wavevector k is 

^k(r) = -—exp(zk.r), (6) 

where V is the volume over which the periodic boundary conditions are applied on the 
wavefunction. Assuming the volume to be a large cube of side L, the wavevectors k take 
the values 

k= ^(riix + n 2 y + n 3 z) (7) 

where ni = 1, 2, 3 . . . with the maximum value such that the magnitude of the largest k is 
kp. The density of k-points in the k-space is therefore ^ and the density of states including 
the spin of the electrons is j-g . Equating the total number of electrons N within the volume 
V to the number of states within a sphere of radius kp, referred to as the Fermi sphere, 
leads to 

k F = (3ir 2 p)^ (8) 

where p = y is the number density of the homogeneous electron gas. Similarly the total 
kinetic energy is calculated by summing the kinetic energy of a state specified by the 
wavevector k over the Fermi sphere. It gives the kinetic energy density or the kinetic energy 
t s per unit volume to be 

ft. 77 



ts = 

IOtt 



2 



= ^(3vr 2 )M (9) 

The local density approximation to the kinetic energy / f (<f>i\ — ^ V 2 \4>%) for the ground- 
state of an inhomogeneous electron gas, such as that in an atom or a molecule, of space- 
dependent density p(r) = Yli fi l<M r )| 2 corresponds to approximating the kinetic energy 
density at each point by the formula above and integrating it over the entire volume. This 
leads to the Thomas-Fermi kinetic energy functional 



T i 0) W = 1^3/ k ^ dY 



= ^(Svr 2 )! J pf(r)dr (10) 

where kp{r) and p(r) at each point in space are related by equation [HJ In equation [TU1 the 
superscript (0) indicates that this is the zeroth-order approximation to the exact kinetic 
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energy for an inhomogeneous electron gas. It is well known to underestimate the exact 
kinetic energy. If the number of up and down spin electrons is different, the functional given 
above can be written in terms of the spin densities p f and p j as 

T[p],pi] = 1 -(T^[2p^] + T^[2pi]) (11) 

Exact kinetic energy for some closed-shell hydrogen-like atoms and the Thomas Fermi ap- 
proximation to it for the same atoms is given in Table I. Table II gives the exact kinetic 
energy for atoms from H to Ne for the density obtained by solving the Kohn-Sham equation 
for it within the Gunnarsson-Lundquist parametrization [?| of the LDA for the exchange- 
correlation energy. As is evident from the Tables, Thomas-Fermi functional underestimates 
the exact kinetic energy by about 5% to 10%. 

The first correction to the Thomas-Fermi functional in terms of the density gradient is 
proportional to the square of the gradient of the density and is given as 3|, |4j 

This is easily derived [l|, [sjj] from the expansion of the response function of a non-interacting 
electron gas. The correction term is also generalized in terms of the spin densities as given 
by equation [Til Equation [12] represents the lowest-order gradient correction to the Thomas- 
Fermi functional. The gradient-corrected kinetic energy + is also given in Tables I 
and II for the atomic systems given there. It is seen that the inclusion of the second-order 
correction brings the approximate kinetic energy closer to the exact one, with the difference 
being less than 1%. The question that we now address is if kinetic energy functionals 
can also be written for excited states using ideas employed to generate functionals for the 
ground-state. 



III. LDA AND GEA UP TO THE SECOND ORDER FOR THE NON- 
INTERACTING KINETIC ENERGY OF EXCITED-STATES 

One straightforward choice for the functionals to be employed for excited-states is to 
use the ground-state functionals described in the section above. However, the way k-space 
is occupied to construct the ground-state functionals does not reflect proper occupation 
of orbitals that are occupied in an excited-state. While for the excited-states of a given 
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number of electrons in a homogeneous electron gas we expect orbitals with relatively larger 
magnitude of wavevectors to be occupied, this does not happen if we use the ground-state 
functional to approximate the kinetic energy of excited-states; as such the ground-state 
functional would underestimate the kinetic energy of excited-states by a much larger amount 
than the proper excited-state functional should. This is shown in Table III where we have 
shown the approximate kinetic-energy, calculated using the ground-state LSD and GEA 
functionals of the section above, of some excited-states of a few hydrogen-like atoms and have 
compared them to the exact kinetic energy. It is seen from the Tables that the ground-state 
functionals indeed underestimate the exact kinetic energy of excited states by significantly 
larger amount than they do for the ground-states. In Table IV, the approximate kinetic 
energies are compared to the exact kinetic energies for excited-state densities obtained by 
solving the Kohn-Sham equation with the Gunnarsson-Lundquist parametrization. Here 
also the error in numbers obtained by applying the ground-state functionals is quite large. 

To develop a kinetic energy functional corresponding to excited-states, we have proposed 
22J in the context of exchange energy that the k-space be split in accordance to the occupa- 
tion of orbitals in the excited-state of a system. This is shown in Fig. 1 for an excited-state 
where some lowest lying orbitals - the core orbitals - are occupied, then there are some 
vacant orbitals and then some more orbitals - the shell orbitals - are occupied. According to 
our method of constructing excited-state functionals for such a system, the k-space is also 
occupied correspondingly with orbitals up to k\ being occupied with k\ given by equation 
orbitals between k\ and k 2 being vacant with k 2 given by equation [3] and again orbitals 
from k 2 to k% being occupied with k% given by equation HI Now steps leading to equation [10] 
are taken to derive the kinetic energy for such a system. The corresponding LDA functional 
2""*(o) f or the excited-state is then given as 

T<°\h, k 2 , h) = ^- 2 j (tf(r) + k!(r) - kl(r)) dr (13) 

The spin-density generalization of equation [13] is given by equation [TT] We call the func- 
tional given by equation [13] the modified Thomas-Fermi functional. In the functional 
above the term J (kf^r)) dr represents the kinetic energy of the core orbitals whereas 
1^2 / (fcf(r) - fcf(r)) dr that of the shell orbitals. Thus £\ (c/>f re | - ±V 2 \<\f° re ) is approxi- 
mated by the former while (4>i hell \ ~ f^ 2 \4>i hel1 ) by the latter. 

We now test the functional of equation [13] for excited-states of hydrogen-like and real 
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atoms. Shown in Table III are the approximate kinetic energies calculated using the func- 
tional of equation [13] for excited-states of hydrogen-like atoms. The numbers shown are 
for excited-states in which orbitals up to principal quantum number ri\ are occupied, those 
from n\ + 1 to n<i are vacant and than n-i + 1 to n 3 are again occupied. It is seen that 
the energies calculated with the functional of equation [JS] are better approximation to the 
exact kinetic energy in comparison to the ground-state functional of equation [TUJ Thus 
while the ground-state kinetic energy functional given by equation [10] underestimates 
the excited-state kinetic energy by a substantial amount, the excited-state functional T*^ 
given by equation [13] has the same accuracy for the excited-states as the ground-state func- 
tional does for the ground-state. In Table IV the numbers for approximate kinetic energy as 
obtained by applying the Thomas-Fermi functional and the modified Thomas-Fermi func- 
tional of equation [13] for the density of Kohn-Sham are given. These densities have been 
obtained by solving the Kohn-Sham equation within the local-spin density approximation 
for the exchange-correlation energy . The numbers are compared with the exact kinetic 
energy £\ /j — |V 2 \4>i) with /j representing the number of electrons in the i th orbital , 
obtained from the Kohn-Sham orbitals \<pi). Similar to the case of hydrogen-like atoms, here 
too the functional of equation [13] gives kinetic energies that are better than those obtained 
from the ground-state functional and has similar accuracy as the ground-state functional 
does for the ground-state densities. We thus conclude that the modified Thomas-Fermi func- 
tional of equation [13] is the correct zeroth-order approximation for the kinetic energies of 
excited-states. More importantly, this indicates that the idea of constructing excited-state 
energy functionals by splitting the k-space is a sound one. 

We next discuss the gradient expansion approximation for the excited-states kinetic en- 
ergy. Since the kinetic energy is a sum of kinetic energy of individual orbitals, the second 
order correction to the kinetic energy for an excited-state can also be written exactly in the 
same manner as the zeroth order approximation given by equation [131 Thus the second-order 
gradient correction to the excited-state kinetic energy is given as 

1 72 J pfr;*,) 72 J p(r;k 3 ) 72 J p(r;* 2 ) 

Like in the functional of equation [13] the term f dr gives the gradient correction 

to the core-orbitals kinetic energy while the last two terms give it for the shell orbitals. 
Here p(r; k) = A is the ground-state density corresponding to Fermi wavevector k. The 
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spin-density generalization of equation [14] is given by equation [TTJ In Tables III and IV 
we also show the second-order corrected kinetic energy 

T * (0 ) + T *(2) for exc ited-states of 
hydrogen-like and real atoms, respectively, and compare it to the exact kinetic energies. 
It is again seen that the second order-correction calculated by using equation O leads to 
improved kinetic energies for the excited-states. 

As an extreme test for the functional of equations [13] and [T4] we apply them to excited- 
states where there are no core electrons, i.e. all the electrons have been excited. The 
exact and approximate kinetic energies for such states are given in Table V for the Kohn- 
Sham densities. Comparison of the numbers given shows the following: while the ground- 
state functional of equations [TD] and [T5] underestimate the exact kinetic energy by very 
large amount, the excited-state functionals of equations [13] and dH bring the error down 
significantly. This again points to the soundness of the idea - that of splitting the k-space - 
behind the construction of these functionals. 



IV. GAZQUEZ-ROBLES FUNCTIONAL FOR EXCITED-STATES 

As mentioned in the introduction, there are other forms of the kinetic energy functional 
for the ground-state that are based on considerations other than the LDA and it gradient 
expansion. One of these approaches constructs a functional by combining the von-Weizsacker 
functional 



rpW 



which is exact for one-orbital systems and the Thomas- Fermi (equation [TO]) functional with 
a correction factor 

where N is the number of electrons in the system. Thus the final functional is 

T (o) 9 [ p ] = T w^ + c( N ) T {o)( p ) ( 17 ) 



The constants Ai = 1.314 and A 2 = 0.0021 for spin-compensated case [4|. It is easily 
generalized to the spin dependent case through equation [TTJ In the functional above, the 
von-Weizsacker term gives accurate kinetic energy for the lowest orbital and the contribution 
from the rest of the orbitals is accounted for by the second term. Thus the factor (l — in 
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the second term plays an important role of subtracting from the Thomas-Fermi functional 
the kinetic energy contribution of the lowest orbital, treated exactly by the first term. 

Applying the same arguments that were used to derive equation [TS] and [T^] - that the 
kinetic energy for an excited-state is written as a combination of the ground-state kinetic 
energy functionals corresponding to the wavevectors k\, k 2 an d k 3 - we write the excited-state 
Gazquez-Robles functional as 

T* g (ki, k 2 , h) = T<®*(p(r- h)) + T(°)f (p(r; k 3 )) - T®*(p(r- k 2 )) (18) 

We have also tested the ground-state Gazquez-Robles functional (equation 17) and its 
excited-state generalization (equation 18) for the excited-states studied in Tables IV and 
V. The results are shown in Tables VI and VII. It is evident from the numbers presented 
that with the Gazquez-Robles functional also, our approach leads to an excited-state func- 
tional that estimates the kinetic energy of an excited-state better than its ground-state 
counterpart. We note, however, that unlike the GEA functional the Gazquez-Robles func- 
tional is not uniformly accurate for all the excited-states studied. This could be because the 
parameters of the functional have been optimized using the ground-state kinetic energies 
of atoms within the Hartree-Fock theory. Nonetheless, by applying our approach to two 
kinetic-energy functionals, which are derived by two different methods, we have shown that 
our method leads to improved functionals for excited states. 



V. CONCLUDING REMARKS 



In this work we have tested the idea of constructing the LDA to excited-state energy 
functionals of time-independent density functional theory by splitting the k-space in the 
context of non-interacting kinetic energy functionals. Our results show that the functionals 
obtained by such a method have the same accuracy for the excited-states as the ground-state 
functionals do for the ground-states. Further, we have shown that gradient correction can 
also be made on such functionals. The general nature of our proposal is evident from the fact 
that applying it to a different kinetic-energy functional also leads to an improved functional 
for the excited-states. In the future we would like to derive the gradient correction given 
by equation [14] in a manner similar to that 35[ for the ground-state, i.e. from the response 
function of the excited HEG. Further, it would also be interesting to see if excited-state 
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functional derived here can be used to approximately calculate excited-state energies by 
employing a variational form for the excited-state densities. 
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the anonymous referees for suggesting that we test our method on functionals other than 
the GEA functional. 
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TABLE I: Closed Shell : Exact and Thomas Fermi (Equation [10] of the text ) and Gradient 
corrected kinetic energy T^°> (Equation [T2] of the text ) for closed shell hydrogen like atoms. 

Numbers given are in atomic units. Percentage errors as shown in brackets under each number 



atoms 


rpExact 


r (o) 


T (0) + T (2) 


He{ls 2 ) 


4 


3.672 


4.134 






(8.2) 


(3.4) 


Be([He]2s 2 ) 


20 


17.719 


19.785 






(11.4) 


(1.1) 


Ne{[Be]2p d ) 


200 


188.849 


202.869 






(5.6) 


(1.4) 


Mg([Ne]3s 2 ) 


304 


284.712 


305.978 






(6.3) 


(0.6) 


Ar([Mg]3p 6 ) 


792 


737.963 


790.652 






(6.8) 


(0.2) 
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TABLE II: Exact Kinetic Energy of ground state of some atoms as obtained by solving the Kohn- 
Sham equation with the Gunnarsson-Lundquist parametrization of the LSD for exchange and 
correlation energy. Numbers given are in atomic units. The exact kinetic energy is compared with 
the Thomas- Fermi (Equation I lOf) and gradient corrected functional (Equation [T3l) 



atoms 


KS 


r (0) 


T (0) + T (2) 


i^ls 1 ) 


0.430 


0.390 


0.438 






(9.3) 


(1.8) 


He{ls 2 ) 


2.780 


2.468 


2.777 






(11.2) 


(1.1) 


Li([He]2s 1 ) 


7.269 


6.521 


7.305 






(10.3) 


(0.5) 


Be([He]2s 2 ) 


14.331 


12.860 


14.347 






(10.3) 


(0.1) 


B([Be]2p l ) 


24.201 


21.649 


24.040 






(10.5) 


(0.7) 


C([Be]2p 2 ) 


37.277 


33.476 


36.980 






(10.2) 


(0.8) 


N([Be]2p 3 ) 


53.899 


48.946 


53.778 






(9.2) 


(0.2) 


0([Be]2p 4 ) 


74.223 


67.084 


73.406 






(9.6) 


(1.1) 


F([Be]2p 5 ) 


98.742 


89.450 


97.472 






(9.4) 


(1.3) 


Ne{[Be]2p % ) 


127.794 


116.838 


126.778 






(8.6) 


(0.8) 
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TABLE III: Kinetic energies (in atomic units) of excited states of hydrogen like atoms. Z gives the 
atomic number of the atom and the excited state is such that the orbitals are occupied upto rti, 
vacant from n\ + 1 to n 2 and occupied again from n 2 + 1 to n.3 and the corresponding approximate 
kinetic energies. The latter are calculated by applying ground-state functionals and 
of Equations [10] and [12] and the excited state functionals of Equations [13] and Q3] The corresponding 
errors are given below each number. 



z 


m 


n 2 


"3 


rji(Exact) 


r (o) 


T (0) + T (2) 




T *(0) +T *(2) 


10 


1 


2 


5 


400 


331.315 


345.737 


389.390 


403.006 












(17.17) 


(13.57) 


(2.65) 


(0.75) 


15 


2 


4 


6 


900 


700.795 


737.727 


873.249 


907.667 












(22.13) 


(18.03) 


(2.97) 


(0.85) 


20 


2 


5 


7 


1600 


1177.696 


1249.214 


1553.078 


1620.010 












(26.39) 


(21.92) 


(2.93) 


(1.25) 


20 


2 


5 


8 


2000 


1486.889 


1558.398 


1952.452 


2019.384 












(25.66) 


(22.08) 


(2.38) 


(0.97) 


25 


3 


4 


7 


3750 


3316.238 


3437.872 


3665.147 


3779.984 












fll 57) 


(8.32) 


(2.26) 


(0.80) 


30 


3 


4 


7 


5400 


4773.266 


4960.075 


5275.695 


5456.973 












(11.61) 


(8.15) 


(2.30) 


(1.06) 


30 


5 


8 


10 


6300 


5410.076 


5597.432 


6171.908 


6353.568 












(14.13) 


(11.15) 


(2.03) 


(0.85) 


35 


2 


4 


6 


4900 


3806.116 


4071.146 


4745.025 


5017.298 












(22.32) 


(16.91) 


(3.16) 


(2.39) 


40 


7 


9 


12 


16000 


14534.00 


14904.93 


15748.793 


16139.184 












(9.16) 


(6.84) 


(1.57) 


(0.87) 


45 


3 


4 


9 


16200 


14521.834 


15010.222 


15880.574 


16426.318 












(10.36) 


(7.34) 


(1.97) 


(1.40) 



1(3 



TABLE IV: Exact kinetic energies (in atomic units) of excited states of some atoms as obtained by 
solving the the Kohn-Sham equation with Gunnarsson-Lundquist parametrization of the LSD for 
exchange and correlation energy and the corresponding approximate kinetic energies. The latter 
are calculated by applying ground-state functionals and + of Equations [10] and [12] 
and the excited state functionals of Equations 1131 and 1 141 The corresponding errors are given below 
each number. 



Atom 


rp(Exact) 


T (0) 


T (0) + T (2) 


j>*(0) 




Be(ls 2 2s°2p°3s 2 ) 


13.768 


12.278 


13.768 


12.459 


13.945 






(10.82) 


(0.0) 


(9.51) 


(1.29) 


O(ls 2 2s 2p 6 ) 


73.094 


64.154 


70.068 


67.545 


73.704 






(12.23) 


(4.14) 


(7.59) 


(0.83) 


O(ls 2 2s 2p°3s 2 3p 4 ) 


65.764 


56.967 


63.516 


59.834 


66.291 






(13.38) 


(3.42) 


(9.02) 


(0.80) 


O(ls 2 2s°2p 3s°3p 6 ) 


65.506 


56.344 


62.815 


59.885 


66.313 






(13.99) 


(4.11) 


(8.58) 


(1.23) 


Ne(ls 2 2s°2p (i 3s 2 ) 


124.508 


109.521 


118.891 


116.152 


125.947 






(12.04) 


(4.51) 


(6.71) 


(1.16) 


Ne(ls 2 2s°2p°3s 2 3p 6 ) 


109.241 


93.430 


103.889 


99.675 


109.920 






(14.47) 


(4.90) 


(8.76) 


(0.62) 


Mg(ls 2 2s°2p 6 3s 2 3p 2 ) 


191.942 


169.083 


182.740 


180.095 


194.392 






(11.91) 


(4.79) 


(6.17) 


(1.28) 


Ar(ls 2 2s°2p 6 3s 2 3p 6 As 2 ) 


501.507 


443.200 


474.770 


474.671 


507.648 






(11.63) 


(5.33) 


(5.35) 


(1.22) 
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TABLE V: Exact kinetic energies (in atomic units) of pure excited states of some atoms (i.e. all the 
electrons have been excited) as obtained by solving the the Kohn-Sham equation with Gunnarsson- 
Lundquist parametrization of the LSD for exchange and correlation energy and the corresponding 
approximate kinetic energies. The latter are calculated by applying ground-state functionals 
and T( ) + T^ 2 ) of Equations [10] and [12] and the excited state functionals of Equations [13] and UM 
The corresponding errors are given below each number 



Atom 


r-p(Exact) 


r (o) 


T (0) + T (2) 




^*(0)_|_^*( 2 ) 


He(2s 2 ) 


0.736 


0.181 


0.263 


0.575 


0.595 






(75.41) 


(64.27) 


(21.88) 


(19.16) 


He(2s°2p 2 ) 


0.676 


0.292 


0.315 


0.614 


0.606 






(56.86) 


(53.43) 


(9.18) 


(10.34) 


Be(2s 2 2p 2 ) 


4.815 


2.160 


2.409 


4.066 


4.079 






(55.14) 


(49.97) 


(15.57) 


(15.30) 


Be{2p 4 ) 


4.565 


2.219 


2.370 


4.337 


4.277 






(51.39) 


(48.09) 


(5.01) 


(6.32) 


Be(3s 2 3p 2 ) 


2.253 


0.475 


0.611 


1.935 


1.943 






(78.94) 


(72.89) 


(14.11) 


(13.78) 


0{2s 2 2p 6 ) 


33.286 


20.073 


21.139 


30.953 


30.781 






(39.70) 


(36.49) 


(7.01) 


(7.53) 


0(3s 2 3p 6 ) 


15.655 


4.504 


5.207 


14.106 


14.077 






(71.23) 


(66.74) 


(9.89) 


(10.08) 


Ne(2s 2 2p 6 3s 2 ) 


60.842 


37.673 


39.607 


57.398 


57.256 






(38.08) 


(34.90) 


(5.66) 


(5.89) 


Mg(2s 2 2p 6 3s 2 3p 2 ) 


98.521 


61.973 


65.046 


93.451 


93.380 






(37.10) 


(33.98) 


(5.15) 


(5.22) 


Mg(2p 6 3s 2 3p 4 ) 


88.796 


51.142 


53.748 


89.023 


88.477 






(42.40) 


(39.47) 


(0.25) 


(0.36) 


Ar(2s 2 2p 6 3s 2 3p 6 As 2 ) 


283.517 


184.194 


192.359 


270.299 


270.662 






(35.03) 


(32.15) 


(4.66) 


(4.53) 
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TABLE VI: Exact kinetic energies (in atomic units) of excited states of some atoms as obtained by 
solving the the Kohn-Sham equation with Gunnarsson-Lundquist parametrization of the LSD for 
exchange and correlation energy and the corresponding approximate kinetic energies obtained by 
applying the ground-state Gazquez functional T^ 9 of Equation [T7] and the excited state functional 
T* 9 of Equation [18] . The corresponding errors are given below each number. 



Atom 


/-p(Exact) 




rp*g 


Se(ls 2 2s°2p°3s 2 ) 


13.768 


14.472 


14.282 






(5.11) 


(3.73) 


0(ls 2 2s°2p 6 ) 


73.094 


69.760 


73.151 






(4.56) 


(0.08) 


O(ls 2 2s 2p°3s 2 3p 4 ) 


65.764 


73.210 


70.675 






(11.32) 


(7.45) 


O(ls 2 2s°2p 3s°3p 6 ) 


65.506 


72.756 


69.242 






(11.07) 


(5.70) 


Ne(ls 2 2s°2p 6 3s 2 ) 


124.508 


118.552 


124.722 






(4.78) 


(0.17) 


Ne(ls 2 2s°2p°3s 2 3p 6 ) 


109.241 


123.323 


118.816 






(12.89) 


(8.77) 


Mg(ls 2 2s°2p 6 3s 2 3p 2 ) 


191.942 


182.374 


192.450 






(4.98) 


(0.26) 


Ar(ls 2 2s 2p 6 3s 2 3p 6 4:s 2 ) 


501.507 


480.681 


508.486 






(4.15) 


(1.39) 
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TABLE VII: Exact kinetic energies (in atomic units) of pure excited states of some atoms (i.e. 
all the electrons have been excited) as obtained by solving the the Kohn-Sham equation with 
Gunnarsson-Lundquist parametrization of the LSD for exchange and correlation energy and the 
corresponding approximate kinetic energies obtained by applying the ground-state Gazquez func- 
tional T^ 9 of Equation [T7] and the excited state functional T* 9 of Equation [18]. The corresponding 
errors are given below each number 



Atom 


/-p(Exact) 


r(o)s 


rp*g 


Be{2s 2 2p 2 ) 


4.815 


2.596 


3.560 






(46.08) 


(26.07) 


Be{2p 4 ) 


4.565 


1.721 


3.303 






(62.29) 


(27.65) 


Be(3s 2 3p 2 ) 


2.253 


1.304 


3.133 






(42.12) 


(39.07) 


0(2s 2 2p 6 ) 


33.286 


14.769 


26.284 






(55.63) 


(21.04) 


0(3s 2 3p 6 ) 


15.655 


7.489 


21.615 






(52.16) 


(38.07) 


Ne(2s 2 2p 6 3s 2 ) 


60.842 


29.171 


51.372 






(52.05) 


(15.57) 


Mg(2s 2 2p 6 3s 2 3p 2 ) 


98.521 


49.498 


86.222 






(49.76) 


(12.48) 


Mg(2p 6 3s 2 3p 4 ) 


88.796 


41.533 


81.642 






(53.23) 


(8.06) 


Ar(2s 2 2p 6 3s 2 3p 6 4s 2 ) 


283.517 


155.167 


265.599 






(45.27) 


(6.32) 
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